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We reexamine effects of tlie p — ui meson mixing mediated by nucieon poiarizations on tfie sym- 
^\ ' metry energy in isospin-asymmetric nuclear matter. Taking into account tlie rearrangement term 

o . 

neglected in previous studies by others, we evaluate the p — lu mixing angle in a novel way within 
, the Relativistic Mean-Field Models with and without chiral limits. It is found that the symmetry 

^ ■ energy is significantly softened at high densities contrary to the finding in earlier studies. As the 

' first step of going beyond the lowest-order calculations, we also solve the RPA equation for the 

00 

' p — uj mixing. In this case, it is found that the symmetry energy is not only significantly softened 

by the p — cj mixing at supra-saturation densities, similar to the lowest-order p — lu mixing, but 
interestingly also softened at subsaturation densities. In addition, the softening of the symmetry 
' energy at subsaturation densities can be partly suppressed by the nonlinear self-interaction of the 

(N 
> 

oo 
cn 



(T meson. 

PACS numbers: 21.65.+f, ll.30.Er, ll.30.Qc, 11.30.Rd 



INTRODUCTION 



Many hadronic and partonic approaches have been used in studying the charge symmetry breaking 



\ (CSB) and its effect on few-body and bulk-matter observables, see, e.g., refs. [1, 2] for reviews. The 

OO . 

' CSB effect can be displayed in the charge-conjugate systems such as the proton-proton and neutron- 

, neutron binary systems whose scattering lengths in the ^5*0 state: a„„ and a^p differ by 10% after 

X 

$H ' deducting the electromagnetic interaction. The CSB can also be used to explain the well-known 

Nolen-Schiffer anomaly of light mirror nuclei. Among the hadronic approaches, the CSB effect has 
been incorporated into the many-body theories by employing explicitly charge-dependent nucleon- 
nucleon interactions, e.g., the Bonn[3], Reid93[4] and V18 potentials [5]. These interactions have been 
adjusted to reproduce the free-space nuclcon-nucleon scattering data. However, the CSB effects in 
free-space and/or symmetric nuclear matter are normally very small. For example, the CSB-induced 
effects in symmetric nuclear matter with the charge-dependent Bonn potential were shown to be 
quite small [3]. The results with the charge-dependent Reid93 potential also showed that the CSB 
effect on the equation of state (EOS) even in isospin-asymmetric nuclear matter is negligible [4]. 

To explore the possible CSB mechanisms at the hadronic level, the p'^ ~lu meson mixing (ROM) [6- 
11] and the hadron mass splitting effects [3, 12] have also been studied extensively. Within the meson 
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mixing picture, the two charge neutral vector mesons undertake a transition to each other through a 
baryon-antibaryon polarization. An essential point is the incomplete cancellation between the oppo- 
site contributions from proton and neutron loops to the transition amplitude [8] . In nuclear medium, 
the polarization acquires an overwhelming enhancement through the particlo-holc excitations. Very 
interestingly, the ROM gets significantly amplified in isospin-asymmetric nuclear matter, see, e.g., 
refs. [10, 13, 14], because of the different numbers of neutrons and protons available. In a nutshell, 
the different stacking of protons and neutrons in the Fermi sea becomes a naissant factor to intrigue 
the CSB in isospin-asymmetric nuclear matter due to the medium-induced ROM. Loosely, we may 
call this the in- medium CSB. 

The ROM in asymmetric nuclear matter also results in significant modifications to properties of 
the isovcctor meson p and its couplings with miclcons [9 11, 13, 14]. Since the potential part of the 
nuclear symmetry energy is dictated by the exchange of p mesons, at least within the Relativistic 
Mean-Field (RMF) models, corresponding modifications to the density dependence of the symmetry 
energy are thus expected. While most earlier studies [10, 11, 15-17] have focused on investigating the 
in-medium meson spectra, a few studies [13, 14] have indeed looked at effects of the medium-induced 
ROM on the symmetry energy. It was found that the p—w mixing angle reaches its maximum of about 
45° in isospin-asymmetric nuclear matter. Moreover, the symmetry energy was sharply stiffened at 
supra-normal densities. However, the rearrangement term, which is crucial for the thermodynamic 
consistency in deriving the pressure, was neglected in these RMF studies [13, 14]. We find that the 
rearrangement term actually causes the pressure to exceed the energy density when the mixing angle 
9 is close to 45°. This certainly violates the causality, namely, the work done by the pressure of the 
system should be less than its total energy. In this work, using the RMF Lagrangian constructed 
by us recently [18, 19] including constraints of a mass dropping scenario according to the Brown- 
Rho (BR) scaling [20-24], we calculate the p — ui mixing angle in asymmetric matter in a novel 
way. We then examine effects of the medium-induced ROM on the in-medium masses of the p and 
Lo mesons and the density dependence of the symmetry energy. As the first step of going beyond 
the lowest-order nucleon polarization, we also evaluate the p — lo mixing within the random phase 
approximation (RPA). Results from the two different levels of approximations will be compared with 
each other. Moreover, in order to examine the model dependence we also perform the calculations 
with the nonlinear RMF model. 

The paper is organized as follows. In Section II, we introduce the medium-induced ROM within 
the RMF model and the novel way to obtain the p — lo mixing angle. We then summarize briefly the 
formalisms for the ROM within both the lowest-order nucleon polarization and the RPA. Results 
and discussions are presented in Section III. A summary is given in Section IV. 
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II. MEDIUM-INDUCED p - uj MESON MIXING WITHIN THE RMF MODEL 

The following lagrangian with the chiral limit serves as the starting point [18, 19]: 



where ■0,(7, V^, and Vp are the fields of the nucleon, scalar, vector, and isovector- vector mesons, 
with their masses M,m*,mQ, and ml, respectively. The meson coupling constants and masses 
with asterisks denote the density dependence, given by the BR scaling with well restrained forms 
and parameters [18, 19]. We note that the BR scaling is still a phenomenological ansatz to mimic 
the properties of the partial restoration of the chiral symmetry although QCD-based effective field 
theories and models support the mass dropping scenario. For instance, as an effective QCD field 
theory, the hidden local symmetry theory has been developed to include the p meson in addition 
to the pion in the framework of the chiral perturbative theory by Harada and Yamawaki and it is 
shown that the p meson becomes massless at the chiral limit [23, 24]. Moreover, the controversy of 
the BR scaling is still unsettled by recent experiments [21], also see the discussions in Ref. [18]. Here 
we resort to the mass dropping scenario according to the BR scaling mainly due to its simplicity. 

The term eV^inV^ in the lagrangian that mixes the isoscalar and isovector mesons breaks the 
charge symmetry in the isospin space. In the free space, the parameter e represents the amplitude 
of the fundamental CSB within the ROM picture. In isospin- asymmetric nuclear matter, since the 
medium effect is dynamically induced, it can be conveniently expressed as 



in terms of a dimensionless parameter 7. The symmetry breaking term can be exorcized at the cost 
of introducing a mixing angle with the unitary transformation [9, 25]: 



£ = ^in^d'' -M + gla - G^j^V^ - G;j^T,Vp^P 



(1) 



e = mom*7, 



(2) 



Vui = 60 sin 6 + a; cos = &o cos ^ — w sin 9, 



(3) 



with 



tan 261 = 



2e 



(4) 



The lagrangian thus becomes 



(5) 
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where both and g* couple differently to the proton and neutron [9]: 

g;^ = G: cos 0-G; sin 0, g*^ = cos + G; sin (6) 

g;^ = G; cos + GZ sin 0, g*^ = G; cos - G^ sin (7) 
The w and p masses are given by 



2 2 cos 26* ' 

:1c2| ^2 ^2 

= 2 2^^^' 

where the masses of p and w mesons shift oppositely with 9, required by the transformational unitarity 
in Eq.(3). The energy density, based on Eq.(5), reads, 

^ ^ i9p^Pp + 9LPnf ^ idppPp - QnpPn? 

2m*^ 2m*2 

where the is embedded in both meson masses and their couplings to nucleons. Given the p — ui 
meson mixing by the nucleon polarization that is beyond the mean field, the rearrangement term 
is inevitably induced to keep the Lorentz invariance of the RMF model. This is a general case in 
the density-dependent RMF models where the density dependence is incorporated beyond the mean 
field, e.g., sec [18, 19, 26]. Besides the contribution from the BR scaling, the density-dependent 
mixing angle results in a new source of the rearrangement term in the pressure 

{9pu,pp + gLpnf , {glppp - glppnf 



2m*2 2mf 



-i<v-.„.+l5;^^"'A|. (10) 

Here, an additional ^-related part —d£/d0-d0/dp neglected in Refs. [13, 14] makes the rearrangement 
term Eq more involved than that in Ref. [18], and the final result is obtained numerically. 

A key point is to derive the p — w mixing angle through the nucleon polarization. In principle, 
one can firstly determine the ROM parameter e or 7 directly from the polarization with the relation 
(ji^j^e = nj^'^(g^) according to the lagrangian (1), similar to the case in the free space [8], and then 
the mixing angle can be obtained through Eq.(4). However, the ROM parameters e and 7 obtained 
in this way have a complicated structure in momentum space [8]. In order to obtain momentum- 
independent ROM parameters in the mean-field approximation, an average procedure is necessary. 
As mentioned in the Introduction, the mixing angle obtained with such an averaging procedure in 
Refs. [13, 14] is large enough to cause the violation of causality. Here, we pursue a new approach. We 




first evaluate the total energy density from the lagrangian (1) but with the explicit ROM replaced by 
the polarization diagram in bulk matter, see Fig. 1. Since the polarization diagram is a contribution 
beyond the RMF approach, the polarization in Fig. 1 is not limited only to its temporal component 
considered in Refs. [13, 14]. The ROM angle appearing in the lagrangian (5) is then determined 
by reproducing this total energy density using Eq.(9). Meanwhile, the ROM parameter e or 7 can 
be calculated through Eq.(4). Nevertheless, here we indeed evaluate the ROM parameter e in the 
lagrangian (1) through the p° — w transition amplitude. In isospin-asymmetric matter, the medium 
effect is dynamically induced due to the incomplete cancellation of the proton and neutron loops. 
In the lowest order of the meson-nucleon coupling, this is equivalent to computing the contribution 
of the lowest-order nucleon polarizations to the energy density. To include the high-order effect, we 
need to solve the RPA (Dyson) equation to sum up relevant nucleon loops to all orders. 



(a) (b) 

FIG. 1: Diagrams for (a) the scattering, and (b) the corresponding contribution to the potential energy 
density. 

The energy density can bo written quite generally as [27-29] £ = T -\- V/2, where the T and V 
are the kinetic and potential parts, respectively. The polarization diagram gives rise to a term V' 
adding to the potential energy density V. Considering V is just a few percents compared to the 
vector potential even in pure neutron matter, we neglect its direct correction to the meson fields 
(or, nucleon self-energies) that are solved in the RMF equations. In fact, a full calculation with the 
meson fields including this term brings about the momentum dependence of the nucleon potential 
which should be fulfilled at least in the Hartree-Fock (HF) framework. This is, however, beyond the 
present RMF treatment. In the evaluation of the V', the polarization tensor is wrapped up by the 
nucleon current only in the Fermi sea. This approach is actually widely used in the study of bulk 
matter [27-29]. The V is thus written as 

- "J (27r)3(27r)3 ^ 2£;^, ^ 2E;^, 



(11) 



where p*° = E*^ = \Jv'^ + M*^, q = P2 — Pi and go = E*^ — E*-^ . The polarization tensor is given 

by 

n^^W = -iGlC; J ^Tr[7^G(fc)7.T3G(fc + q)], (12) 



where G{k) is the nucleon propagator, and the trace is over both the isospin (r) and spin (7) matrices. 
Here, the polarization tensor is given in the lowest order, while the high-order corrections will be 
considered later on by solving the RPA (Dyson) equation. As seen in Eq.(ll), only the real part of 
the polarization is needed to evaluate the V'. Usually, the imaginary part concerns the stability of 
the mean field in small oscillations. Indeed, it also comes into play in the real part of the polarization 
by taking into account the high-order corrections in the RPA equation. First of all, we consider the 
case with the lowest-order polarization. By choosing the frame {qo, q, 0, 0) with q"^ = qQ—q^, one can 
define the longitudinal (L) and transverse (T) components [30]: 

Ul = Has - Hoo, Ut = Hn = U22. (13) 

In free space 11^ = IIt' = Il^{q^). Here, the tensor coupling constant of p meson is not included 
because it is in principle nonrenormalizable. In the medium, the polarization tensor is usually 
decomposed into Feynman and density-dependent parts: H = H^ + H^, while the former is just 
the one in free space but with the in-medium nucleon mass. The trace over the isospin makes the 
Feynman parts of the proton and neutron polarizations to cancel each other, leading to [8] 

n^^'"(«^) = 1,G:g; /'*-(! In ^S^4^^. (14) 



27r2 - ^io Mf-x{l-x)q^ 

In the medium, the longitudinal and transverse components of the density-dependent part are not 
equal [30]: 

np-.^ - -9,n^r*r* \- . f "^'^ e(fc^. - |fc|) E*,,' -k^x' 

' - "--"^^ l^^^ J (2.)3 Et, {q^y-4(k. qf 

d\ Q{kF^ - \k\) q^\^ (x^ -l)/2+ {k-q) 



^=,,n ■> V-'^)' ^li {q^f-Kk-qf ' ^^^^ 

where O is the step function, kp^ with i = p,n are respectively the proton and neutron Fermi 
momenta, k ■ q = E^^qQ — |k|gx, x = cosZk ■ q, and s = ±1 for proton and neutron, respectively. 
Note that only the real parts of the polarization are given here. With the chosen frame for g^, the 
polarization can be written as: 

n^l = -Sm^'H^, ii,u = 1,2. (16) 



Substituting Eqs.(14),(16) in (11) and considering the conservation of baryon number, we obtain 



(27r) Jo J-i [ E^Ai 



Ep.iEp^i I {m*o -q^){ml - q^) 



with Xp = cosZpi • p2. With this expression that surrogates the expUcit ROM, we can obtain 
straightforwardly the total energy density from the lagrangian (1) in the RMF approximation. 

In the above, we calculated the additional energy density from the lowest-order polarization of the 
ROM without including the high-order ones. A more elaborate treatment needs to include the ROM 
and the high-order corrections altogether by solving the RPA equation. The RPA equation for the 
polarizations is given as: 

n,.. = n^. + il^,xV''^Ur.. (i8) 

This equation can be decomposed into the longitudinal and transverse parts as 

ilL = Ul+UlVlUl, (19) 
where the lowest-order polarization matrices are given as: 

li llg lip 

ni^o- ng- , (20) 

-prfjp T-rpUl T-rpp 

and the meson propagator matrices are diagonal [31, 32] 

^T={^^^ n)'T^L=l^ q'DJf I , (21) 



with A = l/{q^ - m*2 + ie), and D^^p = l/{q^ - w^fi + ie). Now, the polarizations of the ROM in 
the RPA can be written explicitly as: 

where the dielectric functions cl and ct are the determinants of the matrices (1 — VlHl) and 
(1 — DtIIt), respectively. Their explicit expressions, together with definitions of the polarizations in 
matrices (20), are given in appendix A. Note that the imaginary parts of the polarizations are also 
included in the RPA equation, and thus finally we take the real part of the 11'"^ in the calculation 
of the energy density. Here, except for 11''", the vacuum part is neglected for all other polarizations 
since we work in the RMF approximation. Substituting above expressions into Eq.(17), we may 
evaluate the additional energy density within the RPA. 

III. RESULTS AND DISCUSSIONS 



Our numerical calculations are carried out based on the SLC [19] and the well-known NL3 [33] 
parameter sets. The SLC parameter set was obtained by reproducing the pressure profile constrained 
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at supra-saturation densities by the collective flow data from relativistic heavy-ion reactions [34] and 
the density dependence of the symmetry energy at subsaturation densities obtained from study- 
ing isospin diffusion in heavy-ion reactions at intermediate energies [35-37], besides the saturation 
properties of nuclear matter. In the following, we first carry out numerical calculations within the 
lowest-order ROM picture. This facilitates easy and direct comparisons with earlier studies by oth- 
ers [13, 14]. In this case, the retardation effect is neglected by taking (/o = as in many other 
calculations [27-29]. This makes the calculations much easier. The small difference in proton and 
neutron masses is retained in our calculations. To go beyond the lowest-order ROM, we have also 
re-calculatcd all key quantities using the RPA. Results of these calculations will be compared with 
those using the lowest-order ROM. 

A. Results with the lowest-order ROM 




FIG. 2: (Color online) The mixing angle (the upper panel) and the parameter 7 (the lower panel) as a 
function of density for different isospin asymmetries and npe matter at /3 equilibrium with the CSLC and 
CNL3. 



The p — LO mixing angle is obtained by reproducing the total energy density including the polariza- 
tion diagram as we outlined in the previous section. Fig. 2 depicts the p — u mixing angle 9 and the 
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dimensionless parameter 7 for different isospin asymmetries a = {Pn — pp) / Pb as a function of density 
with the CSLC (sohd lines) and CNL3 (dashed hnes). The model CSLC is simply the model SLC 
plus the in-medium CSB given by the ROM. Since the in-medium CSB is negligible in symmetric 
matter, both the CSLC and SLC have almost the same EOS for symmetric matter. In order to have 
the same E'sym = 31.6 MeV at saturation density po = 0.16/m~^ [36, 37] as the SLC, Gp (i.e., G* at 
p = 0) in the CSLC is readjusted to 2.82. Similarly, the model CNL3 is constructed from the original 
model NL3 plus the in-medium CSB but with Gp = 3.02 to have the same symmetry energy at 
saturation density as given by the NL3. The 9 is negative which is opposite to that in Refs. [13, 14] 
where only the temporal component of the polarization was considered. Our result is consistent 
with that in Ref. [10] where the 6 was obtained from accounting higher-order contributions in the 
Dyson equation. In free space, the 6 goes to zero in our calculation because the observed energy 
density, actually obtained by subtracting the vacuum expectation value, is zero in the free space. 
The mixing angle is negligibly small for a = 0, and increases with a as one expects. For pure neutron 
matter, the magnitude of the 9 reaches the maximum value of about 40° at very low densities and 
then decreases with the increasing density. Since neutron stars are among the most neutron-rich 
and mysterious objects in the Universe, it is interesting to also examine the mixing angle in the 
neutron-proton-electron (npe) matter at fi equilibrium. As shown also in Fig. 2, the mixing angle 
in the npe matter at equilibrium is generally between those for pure neutron matter and a = 0.5. 
The 7, shown in the lower panel of Fig. 2, varies from almost zero in symmetric matter to about 0.2 
in highly isospin-asymmetric matter at low densities, while at high densities it decreases with the 
CSLC and almost saturates with the CNL3, similar to the mixing angle 9. We stress here that in 
all cases considered, the mixing angle is negative with its magnitude well below 45°. This is rather 
different from that obtained in Refs. [13, 14]. We can thus avoid naturally the unphysical result that 
the p-meson mass becomes imaginary when the magnitude of 6 is close to 45° as seen in Eq.(8). 

The nuclear EOS obtained with the SLC differs mainly in the high density region from that with the 
NL3. The latter also has no nonlinear terms for vector mesons and thus guarantees the simplicity 
of the linear transformation in Eq.(3). It is thus interesting to compare their predictions on the 
mixing angles. First of all, it is seen from Fig. 2 that the maximum magnitude of the mixing angle 
with the CNL3 is also well below 45°. Secondly, it is seen that the difference between predictions 
with the CSLC and CNL3 is mainly at high densities. While the CNL3 predicts approximately 
constant and large mixing angles, the values with the CSLC become almost zero at high densities. 
This is understandable since the mixing angle is originated from the nucleon polarization. In the 
CSLC, the vertices from the polarization are dictated by the BR scaling that plays an important 
role in suppressing the medium-induced ROM contribution as seen in the Eq.(17) and Fig. 1. In the 



10 

CNL3, however, such a suppressing factor does not exist. Therefore, there is a large difference in the 
predicted mixing angle. The difference for the parameter 7 between the CSLC and CNL3, shown in 
Fig. 2, can be understood similarly. 

We now turn to the effects of the medium-induced ROM on properties of vector mesons in asym- 
metric matter. As an example, we consider the npe matter at /3 equilibrium. Fig. 3 displays the 
changes of the vector meson masses versus density in the npe matter at [3 equilibrium for the CSLC 
(upper window) and CNL3 (lower window) models. To separate effects due to the medium-induced 
ROM from those due to the mass dropping scenario, we plot separately the p and m meson masses 
scaled by their respective masses in symmetric matter in the main frame, and the masses scaled by 
their respective masses in vacuum in the CSLC case, which are in fact the BR scaling functions for 
the meson masses [18], in the inset of the upper window. It is seen that the relatively large mod- 
ification occurs mainly within the low density region. Similar to the high density behaviors of the 
mixing angle with the CSLC and CNL3, for the same reason the two models predict slightly different 
masses at high densities. With the CNL3, the splitting of the scaled p and lo masses stays almost 
unchanged at high densities, whereas with the CSLC the splitting vanishes quickly. Nevertheless, 
the modifications to the scaled masses are quite small and remain at the level of about 5 percent or 
less for both models. 




FIG. 3: (Color online) The p and uj meson masses in the npe matter at 13 equilibrium scaled by their 
respective masses in symmetric matter. The inset shows their masses scaled by their respective masses in 
vacuum. $(ps) is the scaling function of meson masses [18]. 
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The mass dropping relative to its vacuum value is a manifestation of the partial restoration of 
chiral symmetry, as suggested in many effective QCD theories. It is worthwhile to mention that some 
evidences for the dropping meson masses was found from studying the dilepton spectra observed at 
the CERN SPS [22, 38] in the 90's. Recently, a downward shift of the oj meson mass was observed at 
the KEK [39] and the ELSA-Bonn [40]. At saturation density, the mass ratios {m.*{a)/mi{Q), i — p,uj) 
with the CSLC are predicted in the npe matter at jS equilibrium {a — 0.92) to be 0.846 and 0.9 for 
the p and uj mesons, respectively. At a = 0.2 which is approximately the average isospin asymmetry 
for the ^'^^Pb nucleus, the scaled masses at saturation density are 0.878 and 0.87 for the p and lu 
mesons, respectively. This is equivalent to a splitting of their in-medium peaks by about 17 MeV. 




Pr/Po 

FIG. 4: (Color online) The symmetry energy as a function of density for three cases (see text). 

To examine effects of the medium-induced ROM on the symmetry energy defined as Esym{p) = 
(|^)p/2p, we compare in Fig. 4 results for three cases: the SLC, CSLCO (SLC plus the additional en- 
ergy density with the medium-induced ROM), and the CSLC. Though this additional energy density 
is negligible in symmetric matter, it becomes important in isospin-asymmetric matter. Consequently, 
the symmetry energy reflecting the cost to go away from symmetric matter is significantly modified. 
It is seen that the effect on the symmetry energy is significant at densities around 1 — 2po- At higher 
densities, the modification fades away because the vector meson coupling constants (see Eq.(17)) 
tends to zero at high densities according to the BR scaling [18]. To see more clearly the influence 
of the medium-induced ROM on the density dependence of the symmetry energy, we also compare 
results obtained with the SLC and CSLC. It is interesting to see that a large softening effect is 
observed at high densities. While considerable progress has been made recently in constraining the 
density dependence of the symmetry energy at subsaturation densities using data from heavy-ion 
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reactions [35-37], experimental information about the symmetry energy at supra- normal densities 
just starts to emerge[41, 42]. The effect of the medium-induced ROM observed here is significant. It 
further adds to the importance to determine experimentally the symmetry energy at supra-normal 
densities. Hopefully, heavy-ion reactions induced by high energy radioactive beams will make this 
possible in the near future [43]. 

Now, let's understand the in-medium CSB more clearly. Firstly, the eflfects of the in-medium CSB 
are related to the modification to the symmetry energy. We may decompose the polarization of the 
ROM into two parts: the density-dependent part and the Feynman (or, the vacuum) part. The 
Feynman part changes the sign of the correction to the energy functional as the charge- conjugate 
operator is acted on, and it would result in a very small contribution to the energy term linear in 
isospin asymmetry. This part is responsible for the usual CSB effect within the ROM picture as one 
performs the charge-conjugate operation. Here we do not see significant numerical contributions to 
the symmetry energy. However, the density-dependent part does not change sign for such an opera- 
tion and mostly contribute to the symmetry energy term quadratic in isospin asymmetry. Secondly, 
the meaning of the in-medium CSB may concern the in-medium nucleon-nucleon interactions. In 
fact, if one discusses the in-medium nucleon-nucleon scatterings, the in- medium CSB has the content 
of the medium-induced ROM, and it is amplified dynamically by the incomplete cancelation between 
the proton and neutron loops in isospin- asymmetric matter. 

Next, we discuss the causality issue in asymmetric nuclear matter. We examine the pressure by 
taking into account the rearrangement term. In addition to the contribution from the density de- 
pendence of the meson masses and coupling constants induced by the BR scaling, the rearrangement 
term in the CSLC includes several ^-entangled terms. In the CNL3, however, the rearrangement 
term is only from a single ^-related source. In Refs. [13, 14], a large 6 close to 45° was obtained in 
isospin-asymmetric matter. In this case, unfortunately, the neglected rearrangement term actually 
dominates the pressure. More specifically, it gives the leading difference between the pressure and 
the energy density, i.e., 

^ " ^ ^ 2mf cos^2g ^g-^" - '^-^'^ d-p' (^^^ 

where rriQ and are taken their respective vacuum values in Refs. [13, 14]. It diverges to positive 
infinity for 6 45° and 89 /dp > [13, 14]. This certainly violates the causality limit of p < £. In 
fact, Eq.(23) remains finite for physical cases. The explicit violation of the causality indicates that the 

extremely large mixing angle obtained in Refs. [13, 14] is clearly inappropriate. It might be useful 
to comment here that the rearrangement term itself, necessarily required by the thermodynamic 
consistency, allows us to check the applicability or rationality of the methodology applied under the 
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extreme conditions of density and isospin asymmetry. To be more quantitative, in our treatment 
we examine the causality condition in pure neutron matter by displaying the relation between the 
energy density and the pressure in appendix B. 

B. Results within the RPA 
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FIG. 5: (Color online) The mixing angle as a function of density for different isospin asymmetries. The 
upper panel is for the CSLC-RPA, and the lower is for the CNL3-RPA. For comparison, the results without 
the RPA correction are also exhibited. 



The RPA equation contains both the real and imaginary parts of the polarizations. As we can 
see, the imaginary part goes away completely if the retardation effect is neglected, namely, by taking 
qo = 0. Thus, to keep the imaginary parts of the polarizations in the RPA equation, we have to 
account for the retardation effect. Here, we define the CSLC-RPAO as the model similar to the 
CSLCO but with the ROM calculated in the RPA. The model CSLC-RPA is the same as the CSLC- 
RPAO but with Gp — 2.8 to reproduce the same symmetry energy at saturation density as the SLC. 
Similarly, we can define the models CNL3-RPA0 and CNL3-RPA (with Gp = 3.36) in reference to 
the CNL3 and NL3. 

Fig. 5 shows the mixing angle with the CSLC-RPA (the upper panel) and CNL3-RPA (the lower 
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panel) for different isospin asymmetries. It is seen that the absolute value of the mixing angle with 
the CSLC-RPA reaches its maximum around the saturation density. Interestingly, this maximum 
is very close to the one obtained from the lowest-order ROM. Away from the maximum point, the 
magnitude of the mixing angle decreases compared to that from the CSLC. This indicates that the 
high-order contributions included in the RPA cancel partially the lowest-order ROM. Concretely, 
this partial cancellation occurs mainly between the longitudinal modes with and without the RPA 
correction since both the transverse mode of the lowest-order ROM and its RPA correction are small 
in the whole density region. It is known that the transverse mode of the lowest-order polarization 
is small compared to the longitudinal one. In the SLC model, the small value of the transverse 
RPA correction is mainly due to the following two facts. Firstly, the small transverse mode of the 
lowest-order ROM does not increase appreciably with the density because of the suppression added 
by the BR scaling. Secondly, the transverse eigencondition ct = can not be satisfied accordingly. 
However, the RPA correction in the NL3 model gives rise to a very different feature at high densities. 
As shown in the lower panel of Fig. 5, the mixing angle with the CNL3-RPA starts increasing from 
about twice the normal density and then becomes positive at density p > 2.5po- Similar to the 
case with the CSLC-RPA, the RPA correction to the longitudinal mode reaches its maximum in 
the vicinity of the saturation density and goes steadily without large changes at high densities in 
the CNL3-RPA calculations. Thus, this very different behavior from that obtained with the CSLC- 
RPA can only be due to the distinct property of the RPA correction to the transverse mode of the 
ROM with the CNL3-RPA. Within the CNL3-RPA, the small transverse mode of the lowest-order 
polarization increases with the density and the transverse eigencondition = can be satisfied at 
high densities (p > 2.0po). This then results in a large RPA correction to the transverse mode of 
the ROM, responsible for the sign change of the mixing angle. Similarly, at high densities the large 
difference between the CNL3 and CNL3-RPA results can be explained mainly by using the transverse 
mode property of the RPA polarization in the CNL3-RPA. Next, we examine the mixing angle in the 
low-density region. At very low densities (p < 0.3po), the mixing angle in highly isospin asymmetric 
matter can be shifted to positive values with both the CSLC-RPA and CNL3-RPA models. This 
shift is due to the reduced Pauli blocking on the various intermediate excitations (polarizations) at 
very low densities. 

Fig. 6 depicts the parameter 7 obtained with the CSLC-RPA and CNL3-RPA. Compared to 
the results shown in Fig. 2, the RPA correction shifts both the maximum value and its position 

moderately in the low-density region for both models. At high densities, larger effects of the RPA 
correction start to appear around 2po. It is seen that the difference in 7 between the CSLC-RPA 
and CNL3-RPA calculations also becomes more appreciable at high densities {p > 2po). These 
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FIG. 6: (Color online) The parameter 7 with the CSLC-RPA and CNL3-RPA for various isospin asymmetries. 



observations are consistent with those from studying the mixing angle shown in Fig. 5. 

Next, wc turn to the RPA correction to the symmetry energy. Shown in Fig. 7 is a comparison of 
the symmetry energies obtained with the ROM with or without the RPA correction. We first examine 
in the upper panel effects of the ROM within the RPA on the symmetry energy with the CSLC-RPAO 
and CSLC-RPA calculations. In the SLC model, consistent with the analysis for the mixing angle 
in the above, the RPA correction to the energy density is dominantly from the longitudinal mode. 
Comparing the results obtained with the CSLC-RPAO and SLC, we can see that with the increasing 
density the medium-induced ROM with the RPA affects only slightly the symmetry energy. This 
is consistent with the effects on the mixing angle shown in Fig. 5. Moreover, by comparing the 
symmetry energy with the CSLC-RPAO to that with the CSLCO, see Fig. 4, we see that the RPA 
correction to the lowest-order ROM is small at high densities. Since the RPA correction reaches 
its minimum at saturation density as shown earlier in Fig. 5, the bump is interestingly formed in 
the CSLC-RPAO results. Noticing also that the RPA correction is almost negligible at saturation 
density, we can infer that the medium-induced ROM with the RPA modifies the symmetry energy 
as much as the lowest-order ROM at saturation density. When the symmetry energy at saturation 
density is constrained with the empirical value, it is then considerably softened at high densities. 
This softening is clearly seen in Fig. 7 by comparing the results obtained with the CSLC-RPA and 
SLC. At high densities, there is almost no difference between the symmetry energies obtained with 
the CSLC and CSLC-RPA. Furthermore, at subsaturation densities the symmetry energy with the 
CSLC-RPA is also softened with respect to the SLC due to the high-order correlations included in 
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FIG. 7: (Color online) The symmetry energy as a function of density including the RPA corrections within 
the CSLC (upper panel) and CNL3 (lower panel). Here, the retardation effects are also included in the 
calculation within the CSLC and CNL3. 

the RPA calculations. 

In the lower panel of Fig. 7, we compare results obtained with or without the RPA correction 
based on the NL3 parameter set. The CNL3-RPA0 and CNL3-RPA results display some noticeable 
bumpy structures. The bump formed at the lower density is associated dominantly with the RPA 
correction to the longitudinal mode of the ROM, since the transverse RPA correction is small at low 
densities. Thus, it can depend sensitively on the nonlinear self-interaction of the a meson through 
the unperturbed propagator A, (see Eqs.(22),(A2)). In the nonlinear RMF model (such as the NL3), 
the cr-meson mass in the unperturbed propagator A becomes the effective mass due to its nonlinear 
self-interaction [32]. Once the self-interaction of the a meson is neglected in the A, we find that the 
bump at the lower density is then shifted upwards and centralized around the saturation density, 
similar to the case with the CSLC-RPAO which is free of the nonlinear meson self-interaction. As 
the nonlinear self-interaction of the cr-meson is necessarily included in the A in calculations with 
the NL3 model, the symmetry energy at subsaturation density is expected to be modified differently 
by the RPA compared to calculations with the SLC model. Noticeably, starting about 2po the 
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symmetry energy becomes softer in calculations with the RPA. Consequently, in the high-density 
region, the difference between the NL3 and CNL3-RPA0 results (lower panel) is significantly larger 
than that between the SLC and CSLC-RPAO results (upper panel). This feature is consistent with 
that shown for the mixing angle and the parameter 7 in Fig. 5 and 6, respectively. Again, this can 
be attributed to the different properties of the RPA correction to the transverse mode of the ROM in 
the CSLC-RPAO and CNL3-RPA0 calculations. Moreover, the lowest-order ROM can largely soften 
the symmetry energy as seen clearly by comparing the results obtained with the NL3 and CNL3. 
This softening by the lowest-order ROM is dictated by the longitudinal mode, while the transverse 
mode plays only a minor role. We note that the RPA correction to the longitudinal ROM suppresses 
considerably the softening due to the lowest-order (longitudinal) ROM. Thus, it is worth stressing 
that in the nonlinear RMF model (NL3) the RPA correction to the transverse mode of the ROM plays 
a crucial role in softening the symmetry energy at high densities. Comparing the symmetry energy 
with the CNL3-RPA to that with the NL3, it is seen that the softening is very significant at high 
densities. Interestingly, the symmetry energy at high densities with the CNL3-RPA is comparable 
to that with the CSLC-RPA which features a typically soft symmetry energy at high densities. 

In this work, we have only considered the medium-induced ROM and its impact on the vector 
meson properties and the symmetry energy. The in-medium u — p mixing may also affect the 
symmetry energy. However, its effect is expected to be much smaller than the ROM due to the small 
a — p polarization in the low-momentum region [44]. Nevertheless, a thorough investigation including 
the G — p mixing may be an interesting topic for a future work. It is also worth noting that in solving 
the RPA equation the Feynman parts of the polarizations except for the ROM are neglected. The 
neglected Feynman parts may be useful for studying the vacuum effect on the symmetry energy. 
This of course needs a lot more efforts in renormalized models. Finally, we also mention that the in- 
medium CSB effect on the EOS of symmetric nuclear matter is negligible. This finding is consistent 
with earlier conclusions as we discussed in the introduction. 

IV. SUMMARY 

In summary, using the RMF models with and without the BR scaling at high densities we studied 
the medium-induced ROM and its effects on the in-medium masses of the p and u) mesons and the 
density dependence of the symmetry energy in isospin-asymmetric dense nuclear matter. The mixing 
angle is obtained by reproducing the additional energy density from the u) — p conversion through 
the polarization diagram. We took into account the rearrangement term. Significantly different 
predictions from some earlier studies are made especially at high densities. We found that the 
symmetry energy is significantly softened instead of being stiffened at high densities by the medium- 
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induced ROM. The p — u mass splitting due to the medium-induced ROM in highly neutron-rich 
matter is appreciable mostly at low densities and is small at high densities. While for the RMF model 
with the BR scaling, the p — ui mass splitting further vanishes at high densities. No large drop of 
the p meson mass is observed in asymmetric matter, contrary to some earlier studies by others. The 
medium-induced ROM is also studied with the RPA. The symmetry energy is significantly softened 
by the ROM in the RPA at both supra-saturation and sub-saturation densities. We also find that 
such a softening at sub-saturation densities due to the oj — p mixing can be partly suppressed by the 
nonlinear self-interaction of the <j meson. 
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APPENDIX A: DIELECTRIC FUNCTIONS 

Prior to giving the dielectric functions used in Eq.(22), we define the nucleon polarizations accord- 
ing to the Feynman rules [30, 45]: 

-—^Tr[iG{k)TaiG{k + q)T,], (Al) 
-oo (27r) 

where the vertices Ta = ig^, — «G^7/n, and —iG*T'yn for a = a,CLj, and p, respectively. For the neutral 
p meson, the vertex is —iG*pTs-f^. For the nucleon propagator, we use the expression in Ref. [45]. 
With the general definition of polarizations, one can give the expression of the H^"^ in Eq.(12). The 
elements of the polarization matrices in Eq.(20) can also be written out explicitly according to the 
above definition. 

The nucleon polarization for the vector meson can be appropriately decomposed into longitudinal 
and transverse modes such as in Eq.(13), (19) and (20). In the RPA equation for the nucleon 
polarization for the vector meson, the longitudinal and transverse dielectric functions in Eq.(22) are 
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given as: 



ei = (i-An'^'^)(i-L>^nr)(i-£'pnf)-Li^£»^(i-An^^)(nf )2 

-2AZ)„z),^nrn7nr, (A2) 



and 



er = (1 - £>c.nr)(l - DpUf^P) - D^Dpin'^^f. (A3) 

For simplicity, in the text we do not give the explicit expressions for all other polarizations except for 
the real part of the II'"^ given since they can be found easily in literatures, e.g., see Refs. [30, 31, 46]. 



APPENDIX B: CAUSALITY CONDITION IN PURE NEUTRON MATTER AT HIGH 

DENSITIES 



In order to examine the causality condition in our calculations, we display in Fig. 8 the pressure 
versus energy density for pure neutron matter. Comparing the results obtained using the CSLC and 
SLC, one sees that the in-medium CSB increases the pressure appreciably in both the low and high 
density regions. The same conclusion but to a lesser extent can be drawn from comparing results 
obtained using the CNL3 and NL3. Unfortunately, our numerical calculations, as illustrated in the 
inset for the high density tail, indicate that the causal condition p < £ is broken with the CNL3 for 
p > 7.2po (here, ro = 0.16fm~^). These features displayed in Fig. 8 can be elaborated by virtue of 
the rearrangement term — Sq p in Eq.(lO). The Sq is given by 

So = Soe + Sosij, (Bl) 

where 

dSde ^ d£ d4>i 

with (f)i being the various scaling functions used in the SLC model [18]. For the CSLC calculations, 
since the 9 approaches zero at very high densities, the J^obr is little changed by the ROM. Thus, 
the increase of the pressure at very high densities with the CSLC compared to the SLC is due to the 
term Sq^. In pure neutron matter, the factor d£/d6 in the Eoe can be derived as 

89 2 COS2 20 ^TO*4 m;4^^» + ^^*2 ^,2 )ynu,ynpf>n- I ) 

In the CSLC, the first term dominates at very high densities because the BR scaling renders the 
effective masses of the vector mesons to be close to zero. With d9/dp> and 9^—0 (see Fig. 2), the 
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rearrangement term —'Eoe P is positive, leading to the increase of the pressure at very high densities. 
On the other hand, since the masses of vector mesons are very large in the CNL3 compared to the 
CSLC, the violation of the causality limit at very high densities is dictated by the second term. We 
note here that the first and second terms of Eq.(B2) have opposite signs for negative 6 values. While 
in the CNL3 calculations, unlike the CSLC case, one has 80/ dp < at very high densities. Thus, 
the second term of Eq.(B2) can increase the pressure and further violate the causality limit in the 
CNL3. It is useful to note that the reason for the causal breakdown here is different from the one in 
Refs. [13, 14]. In the latter, the inappropriate inducement produces a very large mixing angle close to 
45° that is much larger than the one obtained here with the CNL3. Here, the violation of the causality 
limit is actually due to an inconsistent treatment. Considering the whole set of coupled equations 
for the propagators and self-energies, the Lorentz covariance is respected only approximately since 
we neglected the ROM loop-diagram correction to the self-energies. Consequently, this may lead to 
a formal violation of the Hugenholtz-Van Hove (HVH) theorem which is a conservation law ensuring 
the thermodynamical consistency. In isospin-asymmetric matter, the HVH theorem can be written 
as 

-^+p^^ = \[Hn{l+a)+n,{l-a)], (B3) 

where /i„ and /ip are the neutron and proton chemical potentials, respectively. Since the nucleon 
chemical potentials are determined by the nucleon effective Fermi energies and self-energies, the 
neglect of the loop-diagram correction to the self-energies in evaluating the total energy density from 
the lagrangian (1) results in a formal violation of the HVH theorem. It is certainly a drawback of the 
CNL3 calculations. Nevertheless, as shown in Fig. 8, this violation is very weak and only visible at 
very high densities. Also, we notice that this situation is not very rare in calculations using models 
going beyond the mean field level. For instance, it is widely acknowledged that approaches such as 
the relativistic and non-relativistic Brueckner theories also unfortunately encounter similar problems 
in fulfilling accurately the HVH theorem [47-49]. 

In principle, a self-consistent treatment at least in the relativistic Hartree-Fock framework is needed 
to resolve the causality violation at very high densities. However, this is beyond the scope of the 
present work. Here, the loop-diagram correction to the self-energy is estimated to be tiny (roughly 
2% in pure neutron matter), and the slight causal violation occurs only in highly isospin-asymmetric 
matter at very high densities. Therefore, the symmetry energy obtained in the vicinity of symmetric 
matter is not expected to be affected appreciably. Compared to the CNL3 calculations, the causality 
preservation with the CSLC is due to the fact that the mixing angle is almost zero at very high 
densities because of the BR scaling used. As shown in Fig. 8, the CSLC curve ends at the critical 
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density for the chiral restoration [18]. 




FIG. 8: (Color online) The pressure versus energy density in pure neutron matter. The inset shows the 
pressure within the small interval. 

The features of Fig. 8 in the intermediate and low density regions can also be understood. In 
the intermediate density region, one has 89 /dp > within both the CNL3 and CSLC calculations. 
Thus the contributions from the first and second terms in Eq.(B2) cancel out to a large extent 
due to their opposite signs. As shown in Fig. 8, in the intermediate density region, the pressure 
is little modified by the ^-relevant rearrangement term. At low densities, however, a very distinct 
feature is exhibited by comparing the difference between the pressures obtained with the CSLC 
and SLC to that between those obtained with the CNL3 and NL3. In the SLC calculations, the 
pressure is increased significantly by the ROM, while it is just moderately increased in the NL3 
calculations. Since their mixing angles are very close in this density region, the modification from 
the 0-relevant rearrangement term Spe is comparable for both the SLC and NL3 calculations. Thus, 
this large difference can only be due to the rearrangement term — Sosi? P appearing only in the SLC 
parameter set. 



[1] G.A. Miller, B.M. K. Nefkens and I. Slaus, Phys. Kept. 194, 1 (1990). 

[2] R. Machleidt and H. Miither, topical review, J. Phys. G 27, 69 (2001). 

[3] R. Machleidt and H. Miither, Phys. Rev. C63, 034005 (2001). 

[4] G. H. Bordbar, Int. J. Mod. Phys. A18, 3629 (2003); Int. J. Theor. Phys. 43, 399 (2004). 

[5] R. B. Wiringa, V. G. Stoks and R. Schiavilla, Phys. Rev. C51, 38 (1995). 

[6] S.A. Coon, M.D. Scadron and P.C. McNamee, Nucl. Phys. A287, 381 (1977); S.A. Coon and R.C. 

Barrett, Phys. Rev. C36, 2189 (1987); S.A. Coon and M.D. Scadron, tbid, C51, 2923 (1995). 

[7] P.C. McNamee, M.D. Scadron and S.A. Coon, Nucl. Phys. A249, 483 (1975). 

[8] J. Piekarewicz and A. G. Williams, Phys. Rev. C41, R2462 (1993). 

[9] M. Kimura, A. Suzuki and H. Tezuka, Phys. Lett. B367, 5 (1996). 



22 



A. K. Dutt-Mazumder, B. Dutta-Roy and A. Kundu, Phys. Lett. B399, 196 (1997). 

A. K. Dutt-Mazumder, R. Hofmann and M. Pospelov, Phys. Rev. C63, 015204 (2000). 

G. Q. Li and R. Machleidt, Phys. Rev. C58, 1393 (1998); ibid, C58, 3153 (1998). 

M. Kimura, A. Suzuki, K. Suzuki and H. Tczuka Phys. Rev. C56, 3070 (1997). 

A. Ishihara, A. Suzuki, M. Kimura and K. Suzuki, Phys. Rev. C61, 034004 (2000). 

Y. Mori and K. Saito, Phys. Lett. B552, 21 (2003). 

R. M. Aguirre and A. L. be Paoh, Phys. Lett. B603, 13 (2004). 

Y. Muto and K. Saito, Phys. Lett. B659, 565 (2008). 

W. Z. Jiang, B. A. Li and L. W. Chen, Phys. Lett. B653, 184 (2007). 

W. Z. Jiang, B. A. Li and L. W. Chen, Phys. Rev. C76, 054314 (2007). 

G. E. Brown and M. Rho, Phys. Rev. Lett. 66 (1991) 2720. 

G. E. Brown, J. W. Hoh, C.-H. Lee and M. Rho, Phys. Rep. 439, 161 (2007). 

G. Q. Li, C. M. Ko and G. E. Brown, Phys. Rev. Lett. 75 (1995) 4007. 
M. Harada and K. Yamawaki, Phys. Rep. 381,1 (2003). 

M. Harada, hep-ph/0703238. 

M. Binger and C.-R. Ji, Phys. Rev. D60, 056005 (1999). 

H. Lenske and C. Fuchs, Phys. Lett. B345, 355(1995); C. Fuchs, H. Lenske and H. Wolter, Phys. Rev. 
C52, 3043 (1995). 

C. J. Horowitz and B. D. Serot, Nucl. Phys. A464, 613 (1987). 

A. Bouyssy, J.F. Mathiot, V. G. Nguyen and S. Marcos, Phys. Rev. C36, 380 (1987). 
R. Brockmann and R. Machleidt, Phys. Rev. C42, 1965 (1990). 

S. A. Chin, Ann. Phys. 108, 301 (1977). 

C. J. Horowitz and K. Wehrberger, Nucl. Phys. A531, 665(1991). 

S. Reddy, M. Prakash, and J. M. Lattimer and J. A. Pons, Phys. Rev. C59, 2888(1999). 

G.A. Lalazissis, J.Konig and P.Ring, Phys. Rev. C55, 540 (1997). 

P. Danielewicz, R. Lacey and W. G. Lynch, Science 298, 1592, (2002). 

M. B. Tsang, et al., Phys. Rev. Lett. 92, 062701 (2004). 

L. W. Chen, C. M. Ko and B. A. Li, Phys. Rev. Lett. 94, 032701 (2005). 

B. A. Li and L. W. Chen, Phys. Rev. C 72 (2005) 064611. 
G. Agakichiev et al., Phys. Rev. Lett. 75, 1272 (1995). 
M. Naruki et al., Phys. Rev. Lett. 96, 092301 (2006). 

D. Trnka et al., Phys. Rev. Lett. 94, 192303 (2005). 

W. Reisdorf et al. (FORI Collaboration), Nucl. Phys. A 781, 459 (2007). 

Z.G. Xiao, B.A. Li, L.W. Chen, G.C. Yong and M. Zhang, Phys. Rev. Lett. 102, 062502 (2009). 

B. A. Li, Phys. Rev. Lett. 85, 4221 (2000); ibid. 88, 192701 (2002). 

P. Roy, A. K. Dutt-Mazumder, S. Sarkar and J. Alam, J. Phys. G35, 065106 (2008). 

B.D. Serot and J.D.Walecka, Adv. Nucl. Phys. 16, 1(1986). 

K. Lim, and C. J. Horowitz, Nucl. Phys. A 501, 729(1989). 

N. L Kassis, J. Phys. G 6, 205 (1980). 

F. de Jong and H. Lenske, Phys. Rev. C 54, 1488 (1996). 

W. Zuo, L Bombaci and U. Lombardo, Phys. Rev. C 60, 024605 (1999). 



